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We compute the influence action for a system perturbatively coupled to a linear scalar fleld acting 
as the environment. Subtleties related to divergences that appear when summing over all the modes 
are made explicit and clarifled. Being closely connected with models used in the literature, we show 
how to completely reconcile the results obtained in the context of stochastic semiclassical gravity 
when using mode decomposition with those obtained by other standard functional techniques. 
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0\ ■ I. INTRODUCTION 

T-H ■ 

Ch ' The Closed Time Path (CTP) functional formalism has been very useful to study the back-reaction effect in the 

, ^ , context of semiclassical gravity (as well as other aspects of field theory) |jl^-0|. When considering the back-reaction 

problem in semiclassical gravity, one is usually interested only in the gravitational field dynamics whereas the quantum 

>--^ ' matter fields are treated as an environment [SJtI . The results obtained when integrating out the environment degrees 

of freedom are closely connected with the influence functional [g , a statistical field theory method which has proved 

, ' very fruitful to reveal the stochastic nature of open quantum systems (for applications to quantum Brownian motion 

h>. . models see Ref. [g[). In fact, it has been pointed out that semiclassical gravity |g|,^ and effective theories in general 

\^ should exhibit dissipation and noise H . To describe the stochastic character of the system dynamics due to the noise 

C^ ,. induced by the environment, Langevin-type equations are required. Thus, Einstein-Langevin equations have been 

^^ ■ used to address the back-reaction problem in the framework of semiclassical gravity 0,ll^,R,O| . 

^§ When dealing with fields, mode decomposition can be a useful calculational tool since it makes the problem closer 

to quantum mechanical systems (free fields are treated as an infinite set of decoupled harmonic oscillators). The 
main advantage of this method is that the noise and dissipation kernels can be obtained in a rather direct way Il2] 
and, in the context of semiclassical gravity, provides a simple connection with the Bogoliubov coefficients (closely 
related to particle creation effects) [0,0. For each mode no renormalization is required, the need for renormalization 



o 
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I , arises when considering an infinite number of degrees of freedom: it is precisely when summing over all the modes 



that one gets ultraviolet divergences. However, the appearance of distributional functions makes this sum rather 
subtle, the presence of such divergences is not always manifestly evident and misleading results may be obtained. 
In the semiclassical gravity context this is particularly important as one may overlook the need for counterterms to 
renormalize the divergences, which will imply the appearance of finite extra terms when addressing the back-reaction 
problem. These drawbacks do not arise in other treatments based on functional methods typical of quantum field 
C^ ' theory which make no use of mode decomposition, where renormalization seems to be more easily handled [O^JT^]. 
The aim of this report is to show how to reconcile the results obtained by means of a mode-decomposition approach 
with the results based on standard field theory techniques for renormalization in curved space-times |13| . In Sect. O we 
introduce the notation and the model that we are going to work with and evaluate the infiuence action perturbatively. 



A concrete example is considered in Sect. Ill, where sum over modes is performed revealing the appearance of 



divergences, and it is shown how they can be handled. In Sect. [V, the previous results are used to consider models 
treated in the literature which use mode decomposition in the context of stochastic semiclassical gravity and show 
how to reconcile these results with those obtained by usual functional methods. 

II. MODE-DECOMPOSED EXPRESSION FOR THE INFLUENCE ACTION 

To make the description as simple as possible we follow ref. 1q| and consider the whole action for a system described 
by the variable x{t) with action S'[a;(t)] and the environment, described by a free field 4>{t,x) in flat space which has 
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been decomposed in a complete set of modes {uk{x)}: (j){t, x) = J^k 1k{t)uk{x)- The free action S[(f>{t, x)] is local and 
at most quadratic in (f){t,x), since the field is linear. If that is also the case for the term S™*'\x{t),4'{t,x)] describing 
the interaction with the system, the action terms for the environment may be written, after performing the spatial 
integrals and using the completeness relation for the modes, as ^j, S[qk{t)] and ^j, S'™*[x(i), qk{t)] respectively, where 
the action for each mode S[qk{t)\ corresponds to that of a harmonic oscillator. The dynamics is therefore equivalent 
to that of a set of decoupled harmonic oscillators interacting separately with the system. 

5[x(t),^(i,f)] = 5[x(i)] + ^5fe(i)] + ^5r*Wi),?fcW], (1) 

k k 

where S'™*[x(i), gfe(i)] = J dtQk{qk{t)) ■ h{x{t)) with Qk and h being some specific functions. The expression for the 
Feynman and Vernon influence functional g in the interaction picture is: 

F[x+,x_] =e'^^^["+'"-l = /(OinlJjT- (e"'S^"*["-(*)lW+ (e'^^"*["-(*)l) |0 in)^ , (2) 

A: 

where T+ and T~ correspond to the time ordering and anti-time ordering prescriptions respectively. To obtain the 
influence action Sif^ we will treat the interaction term S'™*[x(i)] perturbatively. Taking the logarithm of (0) and 
expanding up to second order in 5'|."*, we get: 

Sif[x+,x-] ~ I dt[G+{t)x+{t) + G-{t)x-{t)] + - I dtdt'[G++{t,t')x+{t)x+{t') + G+-{t,t')x+{t)x^{t') 

+G-+{t,t')x^{t)x+{t') + G—(t,t')x-{t)x^{t')], (3) 

where G+{t) = -G_(t) = Efe(Qfc(ft(i))), G±±(t,i') = Ek ^{{T^Qk{qk(t))Qk{qk{t'))) - {Qk{qk{t))){Qk{qk{t')))) 
and G+^{t,t') = G_+{t',t) = J^k -i{{Qk{qk{t))Qk{qk{t'))) ~ {Qk{qk{t))){Qk{qk{t'))))- All the expectation values 
are considered with respect to the asymptotic in vacuum |0 in)i in the interaction picture. Note that we have integrated 
out the environment degrees of freedom and Sif depends only on the system variables. 

It is important to separate the real and imaginary parts of the influence action because, as is well known ppj , the 
imaginary part is related to the noise that the environment induces on the system, whereas the real part gives the 
averaged dynamics of the system. These are: 
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dt {Qk {qk it))) A (i) + i / dtdt'^ (i) Hk (i, t') A {t') 



^Sif[x+,xJ\ = 

^Sif[x+,x.]=Y, \ fdtdt'A{t)Nk{t,t')A{t') 

J, L '' 



(4) 
(5) 



where we have defined S (t) = ft, {x+ {t)) + h {x- (i)) and A{t) = h {x+ {t)) — h (a:_ (i)) and we have introduced 

Hk{t,t') = Ak{t,t') ~ Dk{t,t') = -2Dk{t,t')e{t-t'), (6) 

which has been expressed in two alternative and equivalent ways for further use. Here the kernels Ai^, D^ 
and Nk are defined as follows: Dk{t,t') = —^{[Qk{qk{t)),Qk{qk{t'))]) is the dissipation kernel and Nk{t,t') = 
\ {{Qk {qk{t)) , Qk {qk{t'))}) — {Qk {qk{t))) {Qk{qk{t'))) is the noise kernel. The dissipation and noise kernels, which 
are related by the fluctuation-dissipation theorem, are antisymmetric and symmetric respectively under interchange 
of t and t'. On the other hand, the kernel Ak {t, t') = | sgn(f — t') {[Qk {qk{t)) , Qk {qk{t'))]) is symmetric and, as we 
will see, it is the part that gives rise to divergences. 

III. SUM OVER ALL THE MODES AND NEED FOR RENORMALIZATION 

For concreteness, let us now consider the case Qk {qk{t)) = ^qkit)'^ {g is a perturbative couphng constant) where 
(j){t, x) is a massless real scalar field satisfying the Klein-Gordon equation in Minkowski space-time. We can use 
the following conventions (note that the label for each mode, fc, corresponds in fact to a three dimensional vector): 
qk (t) = dkfk{t) + dlfl{t), G+{t,t') = {qk{t)qk{t')) = fk{t)n{t') and Gf (t,i') = {Tqk{t)qk{t')) = 0(t-t')fk{t)r,{t') + 
9{t' — t)fk{t')f^{t), where aj, and dk are the creation and annihilation operators for each of the modes Uk{x) in 



which the field (j){t,x) has been decomposed. When properly normalized, fk{t) — (27r) 2{2u>k) ^ exp(— iwfei) with 
^k = (fc^)^. Taking all this into account, we will have 



Dk{t,t') 



Gt{t,t'f-Gt{t',tf 



Ak (t, t') - 4 sgn (t - t') G+ {t, tf - G+ {f, tf 



Nk{t,t') ^-\G+it,tr + Gtit',ty 



(7) 
(8) 
(9) 



To perform the sum over all the modes, we note that we may write Ak + iNk — jG^ and Dk + iNk = ^G^ 
Using the integral representations for Gf and G^, Gl{t,t') = -{2m)-^ j°^^duje-"^^*-^'\uj^ -P +ie)-^ , Gl{t,t') : 
/^ dwe-*"(*-*')(5(w2 - p)e{ijj) we obtain 



F u +/\2 



Gi{t,t') 



(27r) 



dk°^-^k«{t~t') 



dio 



uj^-P + ie 




{lo - k^f -P + ie 



dp e-^'^" ('-'') / duS L^ - P) e {lo) S ({lo - Pf - P) 9 {P ~ lo) . 



(10) 



(11) 



When carrying out the sum over modes ( J^k = ^/(Stt)"^^ / d^^^k ) we note that the G^p term will diverge. Thus, 
we use dimensional regularization to perform the integrals and then expand in powers of (n — 4), where n is the 
spacetimc dimension. The usual procedure gives: 



A{t~t')+iN(t-t') 



V 
32^ 



dfc°e-^"° (*-*') 



n-4 2 



1 f{k''f + ie 



Dit-t')+tN{t-t') ^^ 



dp e-'""" ('-'') e{P) 



(12) 
(13) 



where A = J^k ^k, D = J^k ^k and N — J^k -^k- The second integral is finite and thus, we finally have the finite 
parts: Ae«(i-i') = -{V/32Tr'^) J^^dk°e-"''°<^'-*">\n{kyfi)^, D{t - t') = i(F/327r2) /^^ dfcOe-*'=''(*-*') sgn(fcO) and 
N{t ~ t') = (T//327r)(5(i - t'). The divergent part Adiv [t - t') = -(T//167r(n - A))6{t - t') has been separated in such 
a way that the divergences may be absorbed by counterterms in S[x\. In other QFT contexts (e.g. two interacting 
scalar fields) g,^ the finite contribution from the counterterms can be reabsorbed in the renormalized parameters. 
However, as we will see, in semiclassical gravity some logarithmic finite terms which cannot be reabsorbed arise in 
the counterterms. 



IV. STOCHASTIC SEMICLASSICAL GRAVITY 



As an example we consider the back reaction due to the effect of a small mass or a non-conformal coupling of the 
scalar field (j){x) on a flat Robertson- Walker model [g(lO|l^. We have to make the following substitutions: 



^ (~f^] - — > Q, ( Til 

Qk {qk (i)) • h{x{t)) ^ ^MV?^^^ Hv)) = ^0fe {vf {m' + (C - Cc) R {a (t?))) a (ry)^ 



(14) 



where rj is the conformal time, 0(77) the scale factor, R{a{r])) the scalar curvature, m the scalar field mass and 
^ a dimensionless constant. In those previous works where Einstein-Langevin equations were derived using mode 
decomposition, divergences were not dealt with ||,1^. 

Let us now see how special care is needed with the sum of modes. Take for instance the second definition for Hk in 
(^ and note that, using the real part of (|l^), D{rj — ri') — ^j, Dk{r] — r]') may be written as {V/16tt^) PV(1/(77 — 77')). 
In this case, one would be inclined to write H = J2k ^k ~ D{rj, rj') ■ 9{rj — 77') = (y/167r^) PV(l/(?7 — r]')) ■ 9{i] — 77'), 
but this is an ill-defined product of distributions which may give rise to divergences. A possible way to deal with this 
is by using, instead, the first definition in m) and consider A and D separately: 



H 



Y^Hk=Y.Ak-Y.Dk 



(15) 



where the first term in the last member will be ultraviolet divergent whereas the last term is finite. Now the divergence 
can be clearly identified and one may use the proper counterterm in dimensional regularization to cancel it: 



srHv)] 



32^2 (^_ 4) y ^^v 


-gR' 


327r2 (n - 4) J 


1 dT-^x 


36 

n-A 




+36 - In (a//) - 
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(16) 



The second term in this integral, which is finite, will cause the appearance of extra terms when deriving the Einstein- 
Langevin equation. Using now the results of the previous section, we get total agreement with those results reached 
by functional methods which do not use mode separation |13| . 

A very interesting connection between dissipation and fluctuations in the metric and particle creation has been 
revealed by Calzetta and Hu 0. They computed the energy dissipated by the gravitational field per unit volume 
as pd = J dr]dr]'{dH{rj,ri')/dr])Auj'^{r])Aui'^{r]') (for simplicity we have considered that the asymptotic values of the 
scale factor are ai„ = Uout = 1), and showed that it was equal to the energy density of the created particles, 
p created = (27r)^'^ / 4:TTi''^V i'\(3,^\'^ dv , whcrc Pi, is the Bogoliubov coefficient for the modes with frequency v. However, 

formal use of divergent expressions was made in such a derivation. Our treatment shows clearly that the divergent 



part A{ri,r]') of the kernel H{r],r]') decomposed according to (15) gives no contribution since it is symmetric under 



interchange of rj and rj' and hence the derivative will be antisymmetric: 

Pd = J dW^^l^^lA..^ (^) ^^2 (^,) ^ _ I rf^rf^/^^^Ac.2 (^) A^^ (^') . (17) 

This integral is, therefore, manifestly finite and can be computed using the dissipation kernel obtained from eq. (03) 
thus leading to the same result of Ref. |^ without need to deal with divergent expressions. 
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